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Abstract 

Central to a proof of the Micali- Vazirani general graph maximum matching algorithm 
[MV80] is a definition of blossoms from the perspective of minimum length alternating paths 
[Vaz94]. In this paper we give a new, recursive, definition of blossoms which, although equiva- 
lent to the old one, is more elementary and hence yields new structural facts. These facts, and 
other insights, have yielded a simpler proof of the algorithm. With the benefit of hindsight, 
wc also give a simpler exposition of the algorithm. For all practical purposes, this is still the 
most efficient known maximum matching algorithm. 

1 Introduction 

A proof of the Micali- Vazirani general graph maximum matching algorithm [MV80] was given in 
[Vaz94] ; central to this proof was a definition of blossoms from the perspective of minimum length 
alternating paths. In this paper we give a new, recursive, definition of blossoms which, although 
equivalent to the old one, is more elementary and hence yields new structural facts. These facts, 
and other (more modest) insights, have yielded a simpler proof of the algorithm. Additionally, 
with the benefit of hindsight, we also give a simpler exposition of the algorithm. 

The main new message of [Vaz94] and the current paper is that a substantial graph-theoretic 
development is required not only for proving correctness of the algorithm but also for expounding it 
in a precise manner. As described below, there are some fundamental hurdles to finding minimum 
length augmenting paths efficiently in non-bipartite graphs, a basic step of the algorithm, and 
the algorithm exists precisely because of the structural properties found in this graph-theoretic 
expose. Our contention is that the best way to understand the algorithm is to comprehend how 
it uses this structure as "footholds" to home in on a solution^. 

Matching has had a long and distinguished history spanning more than a century. The following 
quote from Lovasz and Plummer's classic book [LP86] is most revealing: 

Matching theory serves as an archetypal example of how a "well-solvable" problem 
can be studied. ... [It] is a central part of graph theory, not only because of its 
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^Indeed, once the reader is sufficiently comfortable with the structure and the algorithm, they will be able to 
appreciate how the latter blends harmoniously into the former. 



applications, but also because it is the source of important ideas developed during the 
rapid growth of combinatorics during the last several decades. 

Interestingly enough, matching has played an equally central role in the development of the theory 
of algorithms - time and again, its study has not only yielded powerful tools that have benefited 
other problems but also quintessential paradigms for the entire field. Examples of the latter include 
the primal-dual paradigm [Kuh55] , the definitions of the classes P [Edm65b] and # P [Val79] , 
and the equivalence of random generation and approximate counting for self-reducible problems 
[JVV86]. Examples of the former include the notion of an augmenting path [Kon31, Ege31], a 
method for determining the defining inequalities of the convex hull of solutions to a combinatorial 
problem [Edm65a] , the canonical paths argument for showing expansion of the underlying graph 
of a Markov chain [JS89], and the Isolating Lemma [MVV87]. And at the interface of algorithms 
and game theory lies another highly influential matching algorithm: the Nobel Prize winning 
stable matching algorithm of Gale and Shapley [GS62]. 

The running time of the MV algorithm is 0(m^/n ■ a{m, n)) on the pointer model (using Tarjan's 
set union algorithm [Tar75]), where a is the inverse Ackerman function, and 0{m^/n) on the 
RAM model (using Gabow and Tarjan's linear time algorithm for a special case of set union 
[GT85]). Small theoretical improvements to the running time, for the case of very dense graphs, 
have been given in recent years: 0(m-ynlog(n^/m)^) [GK04] and 0{n'^) [MS04], where w is 
the best exponent of n for multiplication of two n x n matrices. The former improves on MV 
for m = n'^~°W and the latter for m = a;(n^'^^); additionally, the latter algorithm involves a 
large multiplicative constant in its running time due to the use of fast matrix multiplication. 
Over the years, other algorithms having the same running time as MV have also been obtained 
[Blu90, GT91]; however, the MV algorithm appears to be the simplest among these. The rest of 
the history of matching algorithms is very well documented and will not be repeated here, e.g., 
see [LP86, Vaz94]. 

The MV algorithm finds augmenting paths in phases; in each phase, it finds a maximal set of 
disjoint minimum length augmenting paths w.r.t. the current matching and it augments along all 
paths. It is easy to show that only 0{^/n) such phases suffice for finding a maximum matching 
[HK73, Kar73]. The remaining task is designing an efficient algorithm for a phase. However, unlike 
the bipartite case, in non-bipartite graphs minimum length alternating paths do not possess an 
elementary property, called breadth first search honesty^ in Section 2. Indeed, in the face of this 
debilitating shortcoming, the problem of finding minimum length alternating paths appears to 
be intractable. It is a testament to the remarkable structural properties of matching that despite 
this, such an efficient algorithm is possible. 

1.1 Overview and contributions of this paper 

This paper is fully self-contained. Sections 2 and 3 give some basic structural properties of 
minimum length alternating paths which are needed for defining blossoms. Especially worth 
mentioning is Theorem 6, which leads to the central notion of base of a vertex. This notion cap- 
tures, in a simple manner, the essential aspect of what would have otherwise been an enormously 
complicated picture; see the diverse-looking examples given in Theorem 6 to illustrate the various 
cases that can arise. 

^Intuitively, it states that one needs to find only shortest paths to nearby vertices in order to find shortest paths 
to far away ones - longer paths are never needed. 
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The new definition of blossoms is presented in Section 4. In turn, blossoms, and their nesting 
structure, impose a strict regimen on how minimum length alternating paths traverse through 
them. These facts, an outcome of the new definition of blossoms, are established in in Section 4. 
These facts go to proving Theorem 25 in Section 5 which contains a crucial structural fact that 
makes the algorithm possible for the non-bipartite case. 

The MV algorithm involves two main ideas: the new search procediire called double depth first 
search (DDFS) and the precise synchronization of events. The former is described in the simplified 
setting of a layered, directed graph in Section 6 and the latter is described in Section 7. Both have 
been presented at a much higher level than in [MV80] and [Vaz94] - without resorting to Pidgin 
Algol code and low level notions, e.g., anomaly edges - and intuitive proofs have been given for 
their correctness. Sections 2, 3, 4 and 5 state key definitions which are essential for understanding 
the algorithm. Knowing the statements of theorems will also help; however, their proofs are 
not essential for this purpose. To the readers who are interested in quickly understanding the 
algorithm, we recommend that they dovetail appropriately between Sections 6 and 7, and the rest 
of the paper. 

Section 8 ties up all the facts and completes the proof. Finally, equivalence of the two definitions 
of blossoms is established in Section 9. It is difficult to overemphasize the importance of well- 
chosen examples for understanding this result; indeed, most of the intuition lies in them and we 
have included several^. 

While writing this paper, we have also thoroughly checked [Vaz94]. Although all the important 
claims made in that paper are correct, not all proofs are complete and "bug- free" . This is especially 
true of Theorems 5 and 6 in [Vaz94], which were proven using complicated case analyses. These 
facts are proven in Theorem 22 in a more principled way using the new definition of blossoms. 
Additionally, we hope that the simplification achieved in this paper will make the result more 
accessible to the research community. For these reasons, the current paper is more suitable for 
pedagogical and archival purposes. 

2 The tenacity of vertices and edges 

All definitions are with respect to a matching M in graph G = (V, E) . We will assume that G 
has at least one unmatched vertex. Throughout, Im will denote the length of a minimum length 
augmenting path in G; if G has no augmenting paths, we will assume that Im = oo. 

Definition 1 (Evenlevel and oddlevel of vertices) The cvenlevel (oddlevel) of a vertex v, 
denoted cvenlevel(i;) (oddlevel(f )), is defined to be the length of a minimum even (odd) length 
alternating path from an unmatched vertex to v; moreover, each such path will be called an 
evenlevel(f ) (oddlevel(u)) path. If there is no such path, evenlevel(f ) (oddlevel(u)) is defined to 
be oo. 

In all the figures, matched edges arc drawn dotted, unmatched edges solid, and unmatched vertices 
are drawn with a small circle. In Figure 1, evenlevels and oddlevels of vertices are indicated; 

missing levels are oo. 

In bipartite graphs, a vertex v can have either an even or an odd length alternating path from an 

^Furthermore, they have been drawn in such a way that they easily reveal their structural properties - this 
involves drawing vertices in layers, according to their minlevel. 
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Figure 1: Evenlevels and oddlevels of vertices Figure 2: Vertex v is not BFS-honest on 
are indicated; missing levels are oo. oddlevel(a) and evenlevel(c) paths. 

unmatched vertex /, not both. Furthermore, minimum length alternating paths are breadth first 
search honest in the sense that if p is a minimum alternating path from f to v and u lies on this 
path then p[f to -u]'^ is a minimum alternating path from f to u^. This elementary property does 
not hold in non-bipartite graphs, e.g., if v lies on the minimum alternating path from / to u, but 
with the opposite parity. This is precisely the reason that the task of finding minimum length 
augmenting paths is considerably more difficult in non-bipartite graphs than in bipartite graphs. 

In Figure 2, v is not BFS-honest on oddlevel(a), oddlevel(5) and evenlevel(c) paths; it occurs 
at length 9 on the first and at length 11 on the other two, even though oddlevel(t;) = 7. Thus 
shortest paths to a, b and c, of appropriate parity, involve longer and longer odd length paths 
to V. Furthermore, this is not just an academic exercise: Suppose this graph had another edge 
(c, d), where d is a new unmatched vertex. Then the only augmenting path uses the evenlevel(c) 
path. Whereas short paths are easy to find in a graph, finding long paths is intractable, e.g. 
Hamiltonian path. Hence, at first sight, the problem of finding minimum length augmenting 
paths, which may involve paths to certain vertices that are much longer than the shortest path, 
appears to be intractable in general graphs. As stated in the Introduction, despite this, the 
remarkable combinatorial structure of matching allows for a very efficient algorithm. 

Definition 2 (Tenacity of vertices and edges) Define the tenacity of vertex v, tenacity(v) = 
evenlevel('u)-|-oddlevel(t'). If {u, v) is an unmatched edge, its tenacity, tenacity(u, v) = evenlevel('u)-|- 

■^This denotes the part of p from / to u. Similarly p[f to u) denotes the part of p from / to the vertex just 
before u, etc. 

^Consequently, an alternating BFS suffices for executing a phase in bipartite graphs, e.g., see Section 2.1 in 
[Vaz94]. 
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evenlevel(f) + 1, and if it is matched, tcnacity(n, t;) = oddlcvcl(n) + oddlevel(w) + 1. 

In Figure 1, the tenacity of each edge in the 5-cyclc is 9 and the tenacity of the rest of the edges 
is oo. Figures 3 and 4 give the tenacity of vertices and edges, respectively, in a more interesting 
graph. 




Figure 3: The tenacity of vertices is indicated; Figure 4: The tenacity of each edge is indicated; 
here a = 13, k = 15, r = 17 and a; = oo. here a = 13, k = 15 and r = 17. 

Lemma 3 The following hold: 

1. If {u,v) is a matched edge, tenacity(u) = tenacity(f) = tcnacity(u, i;). 

2. If {u,v) is an unmatched edge, tenacity(i;) < tenacity(u, i;). 



Proof : If (u, v) is a matched edge, evenlevel(i;) = oddlevel(n) + 1 and evenlevel(ti) = 
oddlevel(f) + 1. And if {u,v) is an unmatched edge, oddlevel(?;) < evenlevel(u) + 1. In each 
case, the lemma follows. □ 

As a result of the first claim of Lemma 3, in several proofs below it will suffice to restrict attention 
to only one of the end points of a matched edge. 

Definition 4 (Limited BFS-honesty) Let p be an evenlevel(i;) or oddlevel(t') path starting 
at unmatched vertex / and let u lie on p. Then \p[f to u] \ will denote the length of this path 

from / to n, and if it is even (odd) we will say that u is even (odd) w.r.t. p. We will say that u 
is BFS-honest w.r.t. p if \p[f to n]| = evenlevel(u) (oddlevel(u)) if u is even (odd) w.r.t. p. 

The usefulness of tenacity is established in Theorem 5, which shows limited BFS-honesty even in 
the non-bipartite case. 
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Observe that in the graph of Figures 2 and 3, tlie vertices ra, b and c arc BFS-lionest on all 
evenlevel and oddlevel paths to the vertices of tenacity a. (However, the vertices of tenacity a 
are not BFS-honest on oddlevel(a) and oddlevel(6) paths.) 

Theorem 5 Let p be an evenlevel(i;) or oddlevel(?;) path and let vertex u £ p with tenacity(n) > 
tenacity(t;). Then u is BFS-honest w.r.t. p. Furthermore, i/ tenacity('u) > tenacity(i;) then 
\p[f to u\\ = minlevel(M). 

Proof : Assume w.l.o.g that p is an evenlevel(v) path and that u is even w.r.t. p (by Lemma 
3). Suppose u is not BFS-honest w.r.t. p, and let q be an evenlevel('u) path, i.e., |g| < \p[f to u\\. 

First consider the case that evenlevel(i;) = maxlevel(u), and let r be a minlevel(i;) path. Consider 
the first vertex of r that lies on p[f to v\. If this vertex is even w.r.t. p then oddlevel(n) < 
|r| + \p[u to u]]. Additionally, cvcnlcvcl(?i) < \p[f to u]\, hence tcnacity(n) < tenacity(f), leading 
to a contradiction. On the other hand, if this vertex is odd w.r.t. p then minlevel(w) = |r| > 
evenlevel(u), because otherwise there is a shorter even path from f to v than evenlevel(t;). We 
combine the remaining argument along with the case that evenlevel(v) = minlevel('i;) below. 

Consider the first vertex, say w, of q that lies on p{u to v] - there must be such a vertex because 
otherwise there is a shorter even path from f to v than evenlevel(f ). If w is odd w.r.t. p then 
we get an even path to v that is shorter than evenlevel(f ). Hence w must be even w.r.t. p. 
Then, q[f to w] op[w to u] is an odd path to u with length less than evenlevel(t;). Again we get 
tenacity(n) < tenacity(i;), leading to a contradiction. □ 

3 The base of a vertex 

Let V be vertex of tenacity t < Im and p be an evcnlcvcl(t;) or oddlcvcl(?;) path starting at 
unmatched vertex /. Define the base of v w.r.t. p, denoted b{p,v), to be the vertex of tenacity 
> t that is furtherest away from / on p. The main fact we will prove in this section is: 

Theorem 6 Let v be vertex of tenacity t < Im- Then the set 

B = {b{p, v) \ p is an evenlevel(i') or oddlevel(?;) path} 

is a singleton. 

We will need some definitions to prove a preliminary fact first. Let p be an alternating path 
starting at unmatched vertex /. Given two vertices u and w on p, if u is further away from / 
on p than w then we will say that u is higher than w. An alternating path that starts and ends 
at matched vertices of p and shares no other vertices of p is said to be a segment on p; clearly, a 
segment will be of odd length. Let s be a segment on p whose endpoints are vertices u and v on 
p. Then we will say that all vertices on p[u to v] are covered by s. A vertex v on p is covered by 
a set S of segments if it is covered by one of the segments in S. A set S of segments on p is said 
to form a flower on p if this set satisfies the following recursive conditions: 

1. S consists of a single segment, s, which starts and ends at even vertices of p. 
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2. S* = S" U {.s}, where S' forms a flower on p and one endpoint of segment s is covered by S' 
and the other is even w.r.t. p. 

3. S" = S" U S"' U {s}, where S' and S" form flowers on p and one endpoint of segment s is 
covered by S' and the other by S". 

In this section, we will denote the matched neighbor of matched vertex u by u'; furthermore, if 
the edge {u, u') lies on p, we will assume that u is even w.r.t. p and u' is odd w.r.t. p. Let S be a 
flower on p and let u and w be the lowest and highest vertices of p covered by S. Then u and w 
will be called the base and tip of the flower, respectively. Observe that u and w will both be even 
w.r.t. p. The length of this flower will be the sum of lengths of all segments in S plus \p[u to 
The proof of the next lemma follows via an easy induction based on the above-stated recursive 
definition of flower: 

Lemma 7 Let S form a flower on p with base u, and let w be any vertex covered by S. Then 
there are odd and even alternating paths from u to w, each of length at most the length of the 
flower. 

Let q be an alternating path that intersects p, possibly at several places. Each odd length 
alternating subpath of q that starts and ends at p will be called a segment of q on p. The proof 
of the following lemma is straightforward. 

Lemma 8 Let q be an alternating path that intersects p such that its first segment starts at even 
vertex w and last segment ends at even vertex u, where w is higher than x. Then w is covered by 
a flower formed by the segments of q on p. 

Let u be a vertex of tenacity t < Im and let p and q be evenlevel(v) and oddlevel(i;) paths, 
respectively, starting at unmatched vertex /. Consider vertices of tenacity > t that appear on 
both p and q; by Lemma 5, each such vertex must be BFS-honest w.r.t. both p and q. Since 
t < ^m, f is one such vertex^. Among these vertices, let b be the highest. We remark that it is 
straightforward to show that neither p(b to v] nor q{b to v] can intersect p[f to b] — otherwise 
either b will have tenacity at most t or there will be a shorter path to v. 

A matched edge {w,w') that lies on both p{b to v] and q(b to v] is said to be a common edge. If 
both these paths traverse this edge in the same direction then we will say that {w, w') is a forward 
edge and otherwise it is a backward edge. A common edge {w, w') is said to be a separator if the 
graph consisting of the vertices and edges of p[6 to t;] U to v] gets disconnected by the removal 
of this edge; clearly, such an edge must be a forward edge. 

We next define three types of frontier edges. Let {w,u/) be a separator edge and let {w,u') and 
{u, u') be the adjacent unmatched and matched edges, respectively, higher up on p. If (w, u') is 
not a separator edge then (t/;, w') is called a separator frontier edge. Next let {w, w') be a matched 
edge that is not a separator edge and such that q{f to w) has no intersection with p{w to v). Then 
if {w,w') is a forward (backward) edge, it will be called a forward frontier (backward frontier). 
Figures 5, 7 and 9, show backward, forward and separator frontiers, respectively, and Figures 6, 8 

®This is precisely the reason for assuming t < Im iri this lemma and beyond; see also the Remark after the proof 
of Theorem 6. 
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and 10 show these three graphs redrawn in a manner that reveals their structure more easily. In 
Figure 9, the evenlevel(u) and oddlevel(t;) paths both use edges (6, w') and {w', w) and the latter 
edge is a separator frontiers. The reason for drawing edge (6, w) is that it provides an oddlevel(u;) 
path. 

The difference between a separator edge and a forward frontier edge is the following: in the 
former case, (6 to v] does not intersect p[f to w) and in the latter case it does. In both cases, 
\p[f to w] \ = \q[f to w]\, since otherwise one of the paths to v can be shortened. 




Figure 5: Edge {w,w') is a backward frontier. Figure 6: 

Lemma 9 Let v be a vertex of tenacity t < Im and let p and q he evenlevel(?;) and oddlevel(?;) 
paths, respectively, starting at unmatched vertex f, such that there are no separators on p{b to v] 
and q{b to v]. Then b{p,v) = b{q,v). 

Proof : We will prove that all vertices on p{b to v] and q{b to v] are of tenacity at most t, 
thereby proving the lemma. By Lemma 3, it suffices to prove this for vertices on p{b to v] that are 
even w.r.t. this path. For each such vertex, say u, we will prove that oddlevel(u) < \q\ + \p[u to v]\ 
thereby proving the claim, since e{u) < \p[f to u\\. 

Let {w,w') be the frontier edge that is closest to {u,u'), such that w is same as or above u. 

First consider the case that {w,w') is a backward frontier. If so, q[f to w] o p[w to u] is an odd 
alternating path which obviously satisfies the length constraint. 

Second consider the case that {w,u/) is a forward frontier. Now, \q[f to w]\ > \p[f to w]\, since 
otherwise we can get a shorter even path to v by splicing the first half of q with the second half of 
p at w. Since {w, w') is not a separator, q[w to v] must traverse an edge of p below {w, w') before 
reaching v. Let {x,x') be the last edge below {w,w') that qlw to v] traverses. Then {x,x') must 
be a backward edge, since otherwise p[f to x] o q[x to v] will be a shorter odd path to v. If (x, x') 
is not the lowest edge of p traversed by q[w to v], then by Lemma 8, x will be covered by a flower 
formed by the segments of q[w to v], and again one can get a shorter odd path to v. Therefore, 
(x, x') must be the lowest edge of p traversed by qlw to v]. By Lemma 8, w is covered by a flower 
formed by the segments of q[w to v]. Let b' be the base of this flower. 
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Now, consider vertex u for which we need to place an upper bound on oddlevel(n). If u hes befow 
a; or if u is not covered by a flower formed by the segments of qlw to v] , then q[f to w] concatenated 
with an odd path from w to b' , through the flower covering w, concatenated with p[b' to u] is an 
odd path to u of the required length. Next assume that u is covered by a flower formed by the 
segments of q[w to v], with the base of this flower being b" . Then p[f to b"] concatenated with 
an odd path from b" to u through this flower is the required odd path; its length is bounded by 

\p\ + \q[w to v]\ < \q\ + \p[u to v]\, 

since to > \p[f to w]\ and we did not use any part of q[f to w]. □ 




Figure 7: Edge {w,w') is a forward frontier. Figure 8: 

Lemma 10 Let v be a vertex of tenacity t < Im and let p and q be evenlevel(u) and oddlevel(t') 
paths, respectively, starting at unmatched vertex f. Then b{p,v) = b{q,v). 

Proof : Lemma 9 dealt with the case that p and q have no separators. Next assume that 
they do and let {w, w') be the lowest separator frontier edge. We claim that w and w' must be 
BFS-honest on p and q - otherwise any evenlevel(w) path must intersect p{w to v] and q{w to v] 
and one can then show that there must be a shorter even or odd path to v than p ov q, respectively. 

Next let us argue that all vertices on p{b to w] and q{b to w] are of tenacity at most t. As in 
Lemma 9, it suffices to restrict attention to vertices on p{b to w] that are even w.r.t. this path. 
For each such vertex, say u, we will prove that oddlevel(M) < \q\ + \p[u to v]\. 

Since w lies on both p and q, and by the definition of vertex b, tenacity(t(;) < t. Let r be an 
oddlevel(u;) path; w.l.o.g. r starts at unmatched vertex /. Since w is BFS-honest w.r.t. p, 
\r\ = tenacity(it;) — \p[f to w]\. Let h be the last vertex of r[f to w) that also lies on p[f to w). 
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First consider the easy case that h = b. If so, r o p[w to u] is an odd path to u of the required 
length. 

Next assume that h ^ b and h is even w.r.t. p. Now, for u on p{h to w), r o p[w to u] is an odd 
path to u of the required length. And for u on p{b to h], q[f to w] o r[w to h] o p[h to u] is the 
required path. Finally, assume that h is odd w.r.t. p. Then, for u on p{b to /i], r op[i/; to u] is 
the required path. And for u on p{h to u;), r op[tf; to tx] is the required path. 

The arguments given above carry over to even vertices on p that lie between any two separator 
frontier edges. The part of p between the highest separator frontier edge and v is similar to the 
no separators case and can be handled using the arguments given in Lemma 9. □ 

Remark: It is easy to construct examples in which {w, w') is a separator edge, and w and w' are 
not BFS-honest on p and q. The claim made in the first para of the proof given above applies 
only to separator frontier edges. 




Figure 9: Edge {w, w') is a separator frontier. Figure 10: 

Proof : (of Theorem 6) In Lemma 10, first fix p and vary q over all oddlevel(u) paths and 
then fix q and vary p over all evenlevel(u) paths to prove that all evenlevel(f ) and oddlevel(f ) 
paths starting at the same unmatched vertex have the same base. Next suppose that there are 
evenlevel('u) and oddlevel('u) paths starting from more than one unmatched vertex. 

Let p and q be evenlevel(u) paths starting from unmatched vertices / and /', respectively. Since 
t < l„i, these paths must meet at a vertex that is odd w.r.t. both paths, say iv' . Now p[w' to v] 
cannot intersect q[f' to u)'), since otherwise we can get a shorter even path to v or an augmenting 
path of length < Z^. Therefore, q[f' to w') op[w' to v\ is also an evenlevel(f) path. Let us name 
it q' . By the fact established in the previous paragraph, b{q',v) = b{q,v). 

Suppose b{p,v) 7^ b{q',v); let b{p,v) = hi and b{q',v) = 62- Clearly, bi lies on p[f to w') and 
62 lies on q'[f' to w'). Therefore, tenacity (it;') < t. It is easy to see that an evenlevel(w') path 
together with p[w' to /] or q'lw' to /'] will yield an augmenting path of length < 1^, leading to a 
contradiction. Therefore, b{p,v) = b{q',v), hence proving the theorem. □ 

Remark: The condition "i < is essential in Theorem 6. To get a counter-example for t = Im, 
consider a graph consisting of an alternating 3-cycle containing unmatched vertex / and a three- 
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length augmenting path from / to /'. In this graph, the tenacity three vertices in the cycle have 
no base. 

Definition 11 (Base of a vertex, and iterated bases) For a vertex v of tenacity t < Im 
define its base to be the singleton vertex in the set defined in Theorem 6. We will denote it as 
base(f ). The iterated bases of v are defined as follows: Define base^(f ) = base(t'), and for A; > 1, 
if tcnacity(base'^(v)) < Im, define base'^"'~^('i;) = base(base'^(i;)). For convenience, we will define 
base°('u) = v even though it is not an iterated base of v. 

Definition 12 (Outer and inner vertices) A vertex v is said to be outer if evenlevel(?;) < 
oddlevel('u) and inner otherwise. 

Clearly, the base of any vertex is outer and hence all iterated bases of a vertex are outer vertices. 
Clearly, basc'^(t;) = v will be inner if v is inner; however, as stated above, base''(?;) is not an 
iterated base of v. In the graph of Figures 2 and 3, the base of each vertex of tenacity a is a, 
tenacity k is b, and tenacity r is c, respectively. Additionally, base^(t;) = a, basc^(i') = b and 
base^(v) = /. Theorems 5 and 6 give: 

Corollary 13 For any vertex v, its iterated bases occur on every evenlevel(?;) and oddlevel(i;) 
path in a BSF-honest manner. 

Definition 14 (Shortest path from an iterated base to a vertex) Let f be a vertex 
of tenacity t < Im and let b be an iterated base of v, i.e., b = base*^(i;) for A; G Z+. By an 
cvenlevcl(5; v) (oddlevel(&; v)) path we mean a minimum even (odd) length alternating path from 
b to V that starts with an unmatched edge. 

Lemma 15 Let v be a vertex of tenacity t < Im and let b = base(?;). Then every evenlevel(i;) 
fbddlevel(v) ) path consists of an evenlevel(6) path concatenated with an evenlevel(6; v) (bddlevel(6; v) ) 
path. 

Proof : Let p be an evenlevel(6) path starting at unmatched vertex / and q be an cvenlevel(6; v) 
path. If their concatenation is longer than evenlcvc^?;) then q must intersect p below b. Let 
{w, w') be the lowest matched edge of p used by q., where w' is even w.r.t. p. Now, using the same 
arguments as those in the proof of Theorem 5, one can show that if w is odd w.r.t. q then there 
is an even path from f to v that is shorter than evenlevel(f ) and if w' is odd w.r.t. q then there 
is a short enough odd path from f to b which gives tenacity(6) <t. □ 

4 Blossoms and their properties vis-a-vis 
minimum length alternating paths 

Definition 16 (Blossoms and their nesting depth) Blossoms will be defined recursively. 
Let b be an outer vertex and t be an odd number such that tenacity (6) > t and t <lm- We will 
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denote the blossom of tenacity t and base b by Bh^t- Define B^^i = and define its nesting depth, 
^QSb,i) = 0. If t > 3, let St = {v\ tenacity(u) = t and base(f ) = 5} and define 

Bj,,t = 5t U j U B,,t-2 

\fe(5tu{6}), V outer 

Define the nesting depth of this blossom, N{Bb^t), to be N{Bb,t-2) if = and N{Bb,t-2) + 1 
otherwise. 

It is obvious from Definition 16 that ii v £ B^^t then tenacity(i;) < t. In the graph of Figures 2 
and 3, the blossom Ba,a consists of vertices of tenacity a, the blossom Bb,K consists of vertices of 
tenacity a and k, and the blossom Bf^r consists of vertices of tenacity a,K and r. The nesting 
depths of these blossoms are 1, 2 and 3, respectively. 

Lemma 17 Let v G Bb,t- Then 3k such that 1 < k < N{Bb,t) md b = base'^(w). Furthermore, 
all the vertices hase{v), . . . , base'^~^('u) belong to Bb,t- 

Proof : By induction on the nesting depth of blossom Bb^t- If N{Bb,t) = Ij by definition, 

b = base(i7). To prove the induction step, suppose N{Bh^t) = ^ + 1. Now, if t; G St, i.e., 
tenacity(t;) = t, then base(t;) = b. Otherwise, 3u G St such that v G Bu.t-2- Clearly N{Bu,t-2) < I 
and base(tt) = b. By the induction hypothesis, 3\k such that I > k > 1 and u = base'^'(t;). Hence, 
b = hase'^'^^{v) and k + 1 < I + 1. Finally, by the induction hypothesis, base(t;), . . . ,hase'^~^{v) 
belong to Bu,t-2- Hence base(ti), . . . , base'^('u) belong to Bb,t- D 

Lemma 18 Let t < t' < tenacity (6) and t' <lm, and let Bb,t and Bb,t' be two blossoms with the 
same base b. Then Bb^t Q ^h,t' ■ 

Proof : 1. The proof is by induction on t' — t. The base case, i.e., t = t' , is obvious. Assume 
the induction hypothesis that Bb^t ^ ^b,t'-2- Now, by Definition 16 it is straightforward that 
Bb,t'-2 C Bb,t'- Hence Bb,t C Bb,t'- □ 

Lemma 19 Let Bb,t be a blossom with base b and tenacity t < Im, and v be a vertex satisfying 
base'^(t;) = b for some k>l. Ift> tenacity(6ase'^~^(v)) then v G Bb,t- 

Proof : The proof is by induction on k. In the base case, i.e., k = 1, base(u) = b. Let 
tenacity(i;) = r, clearly r < t. By Definition 16, v G Bb^r a^id by Lemma 18, Bb^r ^ -^b,*- Hence 

V G Bb,t. 

For the induction step, let base'^~^(?;) = u, tenacity(«) = r < t. By the induction hypothesis, 

V G Bu,r-2- Since base(?x) = 6, by Definition 16, Bu,r-2 ^ Bb^r ^iid by Lemma 18, Bb^r ^ ^b,t- 
Hence v G Bb,t- D 

Lemma 20 Let Bb^t and Bb'^t' be two blossoms such that b G Bb\t'- Then Bb^t C Bb'^t'- 
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Proof : By Lemma 17, there is a > 1 such that V = basc'^'(6) and basc"'^(6), . . . jbase'^ ^(6) G 
By^t'- Clearly, i' > tenacity(base''"^(6)). To prove the statement, we will apply induction on k. 

For the base case, i.e., k = 1, let tenacity(6) = r. Clearly, t < r < t' and I3b,t C By 
Definition 16, C Bf,',r; where the containment is proper since b is not in the first blossom 

but it is in the second one. By Lemma 18, Bb'^r ^ Bjj'^t'- Hence B^^t C Bb',t'- 

For the induction step assume base*^^^(6) = b' . Let base*^(6) = v, and let tenacity(v) = r. Since 

V G Bi)i^t'i ''' ^ t' ■ Clearly, tenacity(basc'^'^"'^(i;)) < r — 2. Therefore, by Lemma 19, 6 G Bv,r-2- 
Furthermore, since base'^(6) = f, by the induction hypothesis, Bh^t C Bv,r-2- 

Since base(i;) = b', by Definition 16, Bv,r-2 Q Bb',r- Since r < t' , Bb',r ^ i3b\t'- Hence, Bh,t C Bb',t'- 
□ 

Theorem 21 The set of blossoms in G forms a laminar family, i.e., two blossoms are either 
disjoint or one is contained in the other. 

Proof : Suppose v lies in blossoms and By^fi. If 6 = 6', we are done by the first claim in 

Lemma 18. Next assume that b ^ b' . Then by Lemma 20, b = base'^(i;) and b' = base'(i;); assume 
k < I. By the second claim in Lemma 20, b = base'^(t;) G Bfji^t'- Finally, by the second claim in 
Lemma 18, Bb^t C Bb\t'- ^ 

For the next theorem, we will need the following definition. Let k = min^j^ G | tenacity(base'(w)) > 
t} and let b = base'^(f). Then we will say that b is the first iterated base of v having tenacity at 
least t. 

Theorem 22 Let tenacity(f ) = t< Im, base(i;) = h, and let p be an evenlevel(6; or oddlevel(6; 
path. Then the following hold: 

1. Let vertex u lie on p with u ^ b and let b' be the first iterated base of u having tenacity at 
least t. If b' lies on p then p[b' ton] is an evenlevel(6'; n) or oddlevel(6'; u) path, depending 
on the parity of \p[b' to u]\. 

2. Let u and b' satisfy the above-stated conditions. Then b' lies on p. 

3. All vertices of p are in B^^t U {h}. 

4- Let vertex u lie on p with u ^ b. Then all iterated bases of u, until b, lie on p. 

Proof : By Lemma 3 we may assume w.l.o.g. that p be an evenlevel(6; v) path. We will prove 
the four claims by strong induction on t. For the base case, let t be the smallest tenacity of a 
vertex in G. In this case all vertices on p, other than 6, will be of tenacity t and will have base b. 
Hence all the claims are satisfied. We prove the induction step below. 

(1). If tenacity(n) = t then b' = b and by Theorem 5, u is BFS-honest w.r.t. p. Now the claim 

follows by Corollary 13 and Lemma 15. 

Therefore if u does not satisfy the claim, t{u) < t. Among the vertices that do not satisfy the 
claim let u be the last one on p. By Lemma 3, u will be even w.r.t. p and by the choice of 
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h\ u € By By the induction hypothesis applied to claim (3), every evenlevel(6'; ■u) path is 

contained in Bb'.t-2 U {b'}. 

Let w be any vertex on p{u to v). Either tenacity(w) = t and base(w;) = 6 or by the choice of u, 
w is of tenacity less than t and satisfies the condition of this claim. In the latter case, let a be the 
first iterated base of w having tenacity at least t. Then, p[a to w] is contained in Ba,t-2 U {a}. 
By Theorem 21, blossom Ba,t-2 is disjoint from Bh'^t-2- Therefore, any evenlevel(6'; -u) path does 
not intersect p{u to v). Hence an evenlevel(5') o evenlevel(6'; u) o p[u to v] path is a shorter path 
than p, leading to a contradiction. 

(2) . If t{u) = t then u is BFS-honest w.r.t. p and by Corollary 13, b' lies on p. Therefore if u 
does not satisfy the claim, t{u) < t. Among the vertices that do not satisfy the claim let u be the 

last one on p. By Lemma 3, u will be even w.r.t. p and by the induction hypothesis applied to 
claim (3), the path cvenlevel(6'; u) is contained in By^t-2 U {&'}• Now there are two cases: 

Case 1: The evenlevel(6') path does not intersect p{u to v\. 

In this case, the path evenlevel(n) o p[u to v] is a shorter path than p, leading to a contradiction. 
Case 2: The evenlevel(6') path intersects p{u to v\. 

Let w be the first intersection of the evenlevel(6') path with p{u to v\. Now w must be even w.r.t. 

p, since otherwise we can obtain a shorter even path to v. Furthermore, tenacity (tw) = t since 
the evenlevel(6') path must enter any blossom of tenacity t — 2 through its base. Let q be the 
path obtained by concatenating the part of the evenlevel(6') path from w to b' together with an 
oddlevel(6'; u) path. Since every oddlevel(u) path consists of an evenlevel(6') path concatenated 

with an oddlevel(6'; path and p[w to n] docs not go through 6', p[w to u] must be longer than 
q. Hence by replacing p[w to u] by g in p yields a shorter path than p, leading to a contradiction. 

(3) . By claim (2), if n is on p and u ^ b, then the first iterated base of u having tenacity at 
least t, say b', lies on p, and by claim (1), p[b' to u] is an evenlevel(5'; or oddlevel(6'; n) path, 
depending on the parity of \p[b' to Therefore, one of the following three must hold: u = b or 
tcnacity(?i) = t and base(u) = b or u £ Bb\t-2 where b' = b oi tenacity(6') = t and base(6') = b. 
In all cases, u £ {Bb,t U {&}). The claim follows. 

(4) . If tenacity = the claim is obvious. Otherwise, tenacity(n) < t and by claim (2), h' lies 
on p. Now by claim (1), p\b' to u] is an evenlevel(6'; u) or oddlevel(6'; tt) path, depending on the 
parity of \p\b' to Hence, we are done by Corollary 13 and the induction hypothesis. □ 

Remark: It is interesting to note that perhaps the most elementary claim in Theorem 22 appears 
to be the following subclaim of (4): If vertex u lies on p, with u ^ b^ then base(«) lies on p. 
However, all our attempts at proving this fact first failed. 

Consider the oddlevel(ra) and evenlevel(c) paths in FigTirc 2. Vertex v lies on both and is not BSF- 
honest w.r.t. either path; the iterated bases of v arc a, b and /. Now, base(a) = b and base(c) = /. 
Let p be an oddlevel(fe; a) path and q be an evenlevel(/; c) path. Clearly, p G Bi),ib U {6} and 
q £ Bf^nVJ {f}. The relevant iterated bases of f , namely a and b lie on p and a, b and / lie on q. 
Finally, p[a to v\ is an evenlevel(a; v) path and q\b to v\ is an evenlevel(6; v) path despite the fact 
that V is not BSF-honest w.r.t. both p and q. 

Theorem 22 leads to the startling fact proved in Theorem 23; this fact also underlines the central 
importance of the notion of base of a vertex. We first need some definitions. Let us say that 
vertex v has all iterated bases defined if for some A; G Z+, base^{v) is an unmatched vertex, and 
vertex v is useful if it lies on a minimum length augmenting path. It is easy to see that any useful 
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vertex has all iterated bases defined; however, the reverse may not hold. Clearly, in both cases, 
tenacity(f ) < Im- 

Theorem 23 Let v be a vertex that has all iterated bases defined and let p be an evenlevel(t;) or 
oddlevel(w) path. Then for each vertex u that lies on p, all iterated bases of u lie on p. 

Proof : First, it is easy to see that vertex u has all iterated bases defined. If u lies on 
p(hasc{v) to v], then by claim (4) in Theorem 22 and Corollary 13, the theorem is true for u. 

Assume that base^{v) is an unmatched vertex, where k € Z"*". liu = base^(y), the claim is clearly 
true. Next assume that u lies on 

p{\:>ase^^^{v) to base'(f)], where 1 < Z < A;. Since base'(t;) is BFS-honest w.r.t. p, 

p[base'"'"^(f ) to base'(t;)] is an evenlevel(base''*"^(i;); base'(v)) path, and again we are done by claim 

(4) in Theorem 22 and Corollary 13. □ 

As an illustration of Theorem 23, consider the evenlevel(c) path, say p, in the graph of Figure 2. 
Vertex v lies on p and so do a = base(v), b = base^('y) and / = base^(u). 

5 Bridges and their support 

Definition 24 (Predecessor, prop and bridge) Consider a minlevel(f ) path and let (u, v) 

be the last edge on it; clearly, (u, v) is matched if v is outer and unmatched otherwise. In either 
case, we will say that ti is a predecessor of v and that edge {u, v) is a prop. An edge that is not a 
prop will be defined to be a bridge. 

The three horizontal edges in Figure 2 are bridges and the rest of the edges of this graph are 
props. 

Theorem 25 Let v be a vertex of tenacity, t < Im, and let p be a maxlevel(t;) path. Then there 
exists a unique bridge of tenacity t on p. 

Proof : Let p start at unmatched vertex /. If i < let base(v) = h. By Lemma 15, p 
consists of an evenlevel(6) path concatenated with an evenlevel(6; v) path. Denote the latter by 
q. lit = Im, let q = p. 

By Theorem 5, each vertex u of tenacity t on g is BFS-honest w.r.t. p. Let us partition these 
vertices into two sets: Si {S2) consists of vertices u such that \p[f to u\\ = minlevel(n) (= 
maxlevel(M)). Let = Si^ {5}. Clearly v & S2. Hence and 52 are both non-empty. Let 
a be the vertex in S-^ having the largest minlevel and c be the vertex in ^2 having the smallest 
maxlevel. Now there are two cases. 

Case 1: a and c are adjacent on p and (a, c) is a matched edge. By Lemma 3, tenacity(a) = 
tenacity(c) = tenacity (a, c). Furthermore, for both a and c, their minlevel is their oddlevel, 
therefore (a, c) is not a prop. Hence it is a bridge of tenacity t. 

Case 2: In the remaining case, a and c must both be outer vertices, and in general, there may be 
several vertices of tenacity less than t between a and c on p. By Theorem 22, the first iterated base 
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of these vertices having tenacity at least t must be a or c. Furthermore, by the third statement 
of Theorem 22, the ones having base a must be contiguous on p and so must be the ones having 
base c. Let a' be the last vertex on p whose iterated base is a; if there is no such vertex, let a' = a. 
Similarly, let c' be the first vertex on p whose iterated base is c; if there is no such vertex, let 
d = c. Then (a', c') will be an unmatched edge of p. We will show that it is a bridge of tenacity 
t. 

By Theorem 22, p[a to a'] is an evenlevel(a; a') path and p[c to c'] is an evenlevel(c; c') path. Now, 
tenacity (o', c') = evenlevel(a') + evenlevel(c') + 1. Substituting evenlevel(c) = t — oddlevel(u), 
evenlevel(a') = evenlevel(a) + evenlevel(a; a') and evenlevel(c') = evenlevel(c) + evenlevel(c; c'), 
we get: 

tenacity(a', d) = evenlevel(a) + evenlevel(a; a') + t — oddlevel(u) + evenlevel(c; c') + 1 = t. 

Clearly, a gets its minlevel from its matched neighbor and if a' ^ a, a' gets its minlcvcl from the 
blossom Ba,t-2- Similarly, c gets its minlevel from its matched neighbor and if c' / c, c' gets its 
minlevel from the blossom Bc,t-2- Therefore {a', d) is not a prop. Hence it is a bridge of tenacity 
t. 

Finally, we show that none of the remaining edges on g is a bridge of tenacity t. Consider an edge 
(i, j) on p\b to a], with i below j on p. If tenacity(j) = t then must be a prop and if t{j) < t 
then j lies in a blossom nested in Bb^t and t{i,j) < t. A similar argument holds for the edges on 
p[c to v]. □ 

Definition 26 (The support of a bridge) Let {u,v) be a bridge of tenacity t < Im- Then, 
its support is defined to be {w \ tcnacity(tt;) = t and 3 a maxlevel(ttj) path containing {u,v)}. 

In the graph of Figures 2, 3 and 4, the supports of the bridges of tenacity a, k and r are the set 
of vertices of tenacity a, k and r, respectively. 



6 Double depth first search 

Double depth first search (DDFS) consists of two coordinated depth first searches. The MV 
algorithm executes DDFS on the given graph, G; however, for ease of exposition we will describe 
it in the simplified setting of a directed, layered graph H. 

The layered graph: The vertices of H are partitioned into h + 1 layers, Ih, ■ ■ ■ lo, with If^ being 
the highest and lo the lowest layer. Each directed edge runs from a higher to a lower layer. Layer 
Ifi contains two vertices, r and g, for red and green. The edges of H are such that from each 
vertex there is a path to a vertex in Iq. A vertex v will be called a bottleneck if every path from 
r to Iq and every path from g to Iq contains v. Observe that v may be in layer Iq. Among the 
bottlenecks, if there are any, the one having highest level will be called the highest bottleneck. If 
so, we will denote it by b. Let 14 (£'{,) be the set of all vertices (edges) that lie on all paths from 
r and g to b. If there are no bottlenecks, there must be distinct vertices ro and go in layer lo and 
disjoint paths from r to ro and g to go- If so, let Ep be the set of all edges that lie on all paths 
from r to ro and from g to go. 

The objective of DDFS: The purpose of DDFS is to find the highest bottleneck if one exists, 
and otherwise, to find distinct vertices ro and go in layer lo and disjoint paths from r to ro and g 
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to go- The running time of DDFS needs to be a funetion of the output in tlie following manner: 
in the former case, DDFS needs to run in time 0{\Ef,\) and in the latter case, it needs to run in 
time 0(|£^p|). In the former case, DDFS also needs to partition the vertices in Vf, — {b} into two 
sets R and G, with r e R and g <E G. Furthermore DDFS needs to find two trees, Tr and Tg, 
rooted at r and g, such that the set of vertices visited by them is i?U{6} and GU{b}, respectively. 
The two trees are found by the two DFSs, called red DPS and green DFS, respectively. 

The two DFSs and their coordination: The two DFSs maintain their own stacks, Sr and Sg, 
which start with r and g, respectively. At any point in the search, each stack contains the vertices 
that have been visited by the corresponding DFS but have not yet backtracked from. Each DFS 
performs the usual operations, with one important modification. The latter is described below 
when we give the rules for coordinating the two DFSs. Because edges in H go from higher to lower 
levels, neither DFS will encounter back edges. The two DFSs start with their center of activity 
at r and g, respectively. Assume that the center of activity of a DFS is at u and it searches along 
edge {u,v). If v is not yet visited by either search, it pushes v onto its stack and moves its center 
of activity to v. In this case, u is declared parent of v. If v is already visited by either search, 
it considers the next unsearched edge incident at u - see below an exception, which is also the 
important modification mentioned above. If all edges incident at u have already been searched, 
it pops u from its stack and moves the center of activity to the parent of u. 

We next give the rules for coordinating the two DFSs. To meet the running time requirement, we 
posit that if b is the highest bottleneck in H, then neither DFS will search along any edges below 
b. This gives our first rule: as soon as the center of activities of the two DFSs are at different 
levels, the one that is higher must move to catch up. If both are at the same level, we adopt the 
convention that red moves first. The exception mentioned above happens when one DFS searches 
along edge {u, v) and v happens to be the center of activity of the other DFS. In this case, v could 
potentially be a bottleneck and the two DFSs first need to determine whether it is. Furthermore, 
if V is not a bottleneck, the two DFSs need to determine which tree gets v. 

When the two DFSs meet at a vertex: The procedure followed by the two DFSs at this 

point is the following, independent of which one got to v first. Let us assume that the red and 
greed DFSs reached v via edges {vr,v) and {vg,v), respectively. First the green DFS backtracks 
from V and tries to reach a vertex, say w, with w ^ v and \e\e\{w) < level(t;). If green succeeds, 
red moves its center of activity to v and it adds v to R and edge (ur, v) to T^, and DDFS proceeds. 
If the green fails, its stack, Sg, must be empty. Next, the red DFS backtracks from v and tries 
to reach a vertex, say with w ^ v and \eve\{w) < level(t;). If red succeeds, green moves its 
center of activity to v; however, it does not push v onto Sg, since it has already backtracked from 
V. In addition, it adds v to G and edge {vg,v) to Tg, and DDFS proceeds. If red also fails, then 
its stack also must be empty and v is indeed the required bottleneck. If so, v is added to both 
R and G and edges {vr,v) and {vg,v) are added to T,. and Tg, respectively, and DDFS halts. If 
DDFS does not find a bottleneck, then eventually the two DFSs must find distinct vertices in Iq. 

Theorem 27 DDFS accomplishes the objectives stated above in the required time. 

Proof : The main difference between a usual DFS and the two DFSs implemented in DDFS 

arises when the two DFSs meet at a vertex, say v at layer Ij. Observe that once one DFS reaches 
a vertex at layer I, say, at every future point, the center of activity of at least one DFS will be at 
level I or lower. Therefore, since both DFSs just moved from higher layers to layer Ij, no other 
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vertices at layer Ij or lower have yet been explored. Therefore, if v is not a bottleneck, there is an 
alternative path that reaches layer Ij or lower and which has not been explored so far. Since at 
this point both DFSs will consider all ways of finding such an alternative, they will find one and 
DDFS will proceed. On the other hand, if v is a bottleneck, there is no such alternative, and after 
considering all possibilities, both stacks will become empty and v will be declared the bottleneck. 

If there is no bottleneck in H, by the arguments given above, the two DFSs will not get stuck 
at any vertex. Hence, one of them must reach a vertex at layer Iq, say v. Since v is also not a 
bottleneck, even if the two DFSs meet at v, one of them will find a way of reaching another vertex 
at layer Iq. 

Clearly, the only edges searched by DDFS are those in Ef, in the first case and Ep in the second. 

Furthermore, each such edge is examined at most twice, once in the forward search and once 
during backtrack. Hence DDFS accomplishes the stated objectives in the required time. □ 

7 The algorithm 

The MV algorithm can be viewed as an intertwining of an alternating BSF (similar to the one 
used for bipartite graphs) with a number of DDFSs. The first part of the algorithm for a phase 
finds the evenlevels and oddlevels of vertices and marks the graph in an appropriate manner. 
This part is organized in search levels, the first one being search level 0. Let jm = ilm — 
where Im is the length of a minimum length augmenting path in G. Then during search level 
jm, a maximal set of augmenting paths of length 1^, is found and the current phase terminates. 
If Im = oo, i.e., there are no augmenting paths, the algorithm will reach a search level at which 
it has found the minlevels and maxlevels of all vertices reachable via alternating paths from the 
unmatched vertices. At this point it will halt and output the current matching, which will be 
maximum. 

7.1 Synchronization of events, and procedures MIN and MAX 

As stated in the Introduction, a key idea in the MV algorithm is the precise synchronization of 
events - we describe this next. During search level i procedure MIN finds all vertices, having 
minlevel(?;) = i + 1 and assigns these vertices their minlevels. For each edge MIN scans, it also 
determines whether it is a prop or a bridge. After MIN is done, procedure MAX finds all vertices, 
V, having tenacity('u) = 2i + 1 and assigns these vertices their maxlevels; their minlevels are at 
most i and are already known. See Algorithm 1 for a summary of the main steps. 
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Algorithm 1 At search level i: 

1. MIN: 

For each level i vertex, u, search along appropriate parity edges incident at u. 

For each such edge {u,v), if {u,v) has not been scanned before then 
If minlevel(v) ^ i + 1 then 
minlevel('u) ■<— i + 1 
Insert u in the list of predecessors of v. 
Declare edge {u,v) a prop. 
Else declare {u, v) a bridge. 
End For 

End For 

2. MAX: 

Find the support of each bridge of tenacity 2i + 1 using DDFS. 
For each vertex v in the support: 

maxlevel(v) ■<— 2z + 1 — minlevel(v) 

End For 



In each search level, procedure MIN executes one step of alternating BFS as follows. If i is even 
(odd), it searches from all vertices, u, having an evenlevel (oddlevel) ofi along incident unmatched 
(matched) edges, say (n, v). If edge {u,v) has not been scanned before, MIN will determine if it 
is a prop or a bridge. If v already been assigned a minlevel of at most i, then {u, v) is a bridge. 
Otherwise, v is assigned a minlevel of i + 1, tt is declared a predecessor of v and edge {u,v) is 
declared a prop. Note that if i is odd, v will have only one predecessor - its matched neighbor. 
But if i is even, v may have one or more predecessors. 

The algorithm also constructs for each odd number t, the list of bridges of tenacity t, B{t). For 
each edge that is declared a bridge, as soon as the appropriate levels of its two cndpoints are 
known, as stated in Definition 2, its tenacity is ascertained and it is inserted into the appropriate 
list. Lemma 31 below proves that by the end of execution of procedure MIN at search level i, the 
algorithm would have determined the tenacity of all bridges of tenacity 2i + l. 

At this point, procedure MAX uses DDFS to find the support of each bridge in the set B{2i + 1). 
This yields all vertices of tenacity 2i + 1 and their maxlevel can be calculated as indicated in 
Algorithm 1. Knowing maxlevels of these vertices helps the algorithm determine the tenacity of 
bridges incident at them. 



19 




Figure 11: A new petal- node is created after Figure 12: Vertices u and v are in the support 
DDFS on bridge (ri, r2). of both bridges of tenacity of 11. 

7.2 Running DDFS on G 

Wc now describe how DDFS is executed on the given graph, G; as stated in Section 6, the graph 
if is a considerably simplified version of G. For graph H, we specified its vertices, including 
a layer for each vertex, and its edges. We will specify the vertices and edges gradually below; 
however, at the outset we state that the layer of any vertex will be its minlevel. 

In the graph of Figure 11, MAX will perform DDFS on bridge (ri, r2), of tenacity 9, during search 
level 4, by starting two DFSs at vertices ri and r2, respectively. We first state a preliminary rule 
for determining the edges of the corresponding graph H - the preliminary rule will suffice for this 
first DDFS. If the center of activity of a DFS is at u then it must search along all edges (u, v) 
where v is a predecessor of u. 

Clearly, this DDFS will terminate with the bottleneck h. It will visit the four vertices which 
constitute the support of bridge (ri,r2); observe that h is not in the support of bridge (ri,r2). 
These four vertices form the blossom 65.9. However, in general, DDFS may not find an entire 
blossom but only a part of it. The set of vertices it does find will be called a petal; a blossom, in 
general, is a union of petals. Thus the four vertices of tenacity 9 are said to belong to the new 
petal. 

At this point, the algorithm creates a new node, called a petal-node; this has the shape of a donut 
in Figure 11. The four vertices of the new petal point to the petal- node; to avoid cluttering Figure 
11, only one vertex is pointing to the petal- node. The bottleneck, 6, is called the hud of the petal. 
The new petal-node points to the two endpoints of its bridge, r\ and r2, and to its bud, 6; the 
reason for the former will be clarified in Section 7.3 and for the latter below. 

Definition 28 (The bud of a vertex) If vertex i; is in a petal and the bud of this petal 
is h then will say that the hud of v is b, written as bud(t;) = 6; bud(v) is undefined if v is not 
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Figure 13: DDFS is performed on the left bridge Figure 14: DDFS is performed on the right 
first, then right. bridge first, then left. 

in a petal. Next we define the function bud*(i;). If v is not in any petal, bud*(t;) = v else 

hud*{v) = bud*(bud(t;)). 

Observe that the bud of a petal is always an outer vertex. Unlike basc(t;) which was defined graph- 
theoretically and independent of the run of the algorithm, bud(f ) depends on the manner in which 
the algorithm breaks ties; for examples, see below. Additionally, bud*(v) not only depends on the 
particular run of the algorithm, it keeps changing as the algorithm proceeds. At any point in the 

algorithm, its latest value is used. 

In Figure 11, a second DDFS is performed on bridge (/i, I2), of tenacity 11, during search level 5. 
To describe this DDFS, we need to state the complete rule for defining edges of the corresponding 
graph H: if the center of activity of a DFS is at u and it searches along edge {u,v), where v 
is a predecessor of u, then DFS must move the center of activity to bud*(i;). Thus, when the 
DFS which starts at I2 searches along edge (hj^i), it moves the center of activity to b. It does 
so by following the pointer from ri to its petal-node and from the petal-node to the bud of this 
petal-node. In the process, the center of activity has jumped down more than one layer. The 
edges of H were allowed to jump down an arbitrary number of layers in order to model this. 

This DDFS will end with bottleneck /. The new petal is precisely the support of bridge (hjh) 
and consists of the eight vertices of tenacity 11 in Figure 11, which includes b. Once again, a 
new petal-node is created and these eight vertices point to it. In addition, the petal-node points 
to li,l2 and to /. Observe that the blossom Bf^n consists of these eight vertices and the four 
vertices of blossom B^^g. 

Next consider the graph of Figure 12 which has two bridges of tenacity 11, (^1,^2) and (ri,r2). 
Observe that vertices u and v are in the support of both these bridges. Hence, the support of 
bridges need not be disjoint. Observe also that the base of these vertices is not a but b; note that 
tenacity(a) = 11. Moreover, there is only one blossom in this graph, i.e., Bb^u. 
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MAX will perform DDFS on these two bridges in arbitrary order during search level 5. Figure 
13 shows the result of performing DDFS on (^1,^2) before (ri,r2). The first DDFS will end with 
bottleneck a. The new petal is precisely the support of (^1,^2), consisting of six vertices of tenacity 
11, including u and v. The second DDFS, performed on (ri,r2), will end with bottleneck b and 
the new petal is precisely the difference of supports of (ri, r2) and (li, I2), i.e., the remaining eight 
vertices of tenacity 11, including a. During this DDFS, when the DFS starting at ri searches 
along edge (ri, u), it realizes that u is already in a petal and it jumps to a, the bud of this petal. 
This ensures that a vertex is included in at most one petal. In general, at the end of DDFS 
on bridge {u,v), the new petal will be the support of {u,v) minus the supports of all bridges 
processed thus far in this search level. 

Figure 14 shows the result of performing DDFS on (ri,r2) before (hjh)- The first petal is the 
support of (ri,r2), i.e., 10 vertices of tenacity 11, including a, u and v. The second petal is the 
difference of supports of (Zi, Z2) and (ri, r2), i.e., 4 vertices of tenacity 11. In both cases, the union 
of the petals found is the blossom Hfe^n. 

Figure 15 illustrates the second way in which DDFS may end, i.e., instead of a bottleneck, it finds 
two unmatched vertices; this happens when DDFS is performed on bridge {u,v) of tenacity 11. 
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Figure 15: DDFS on the bridge of tenacity 11 Figure 16: DDFS performed on bridge (n,u) 
ends with the two unmatched vertices. starts the two DFSs at a and 6, respectively. 

We need to point out one final rule: if DDFS is performed on bridge (ti, v), the centers of activity 
of the two DFSs must start at bud*(u) and bTid*(i'). This rule was vacuous so far, but will 
be applicable while processing bridge [u, v) in Figure 16. The tenacity of this bridge is 19 and 
it will be processed by MAX in search level 9. At that point in the algorithm, the bridges 
of tenacity 15 and 13 would already be processed and u and v will already be in petals, with 
bud(u) = bud*(n) = a and bud(v) = bud*(v) = h. Hence the centers of activity of the two DFSs 
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will start at a and b, respectively. Similarly, in Figure 19, when DDFS is performed on bridge 
{u,v) of tenacity 15, bud*(ti) = u and bud*(t;) will be the unmatched vertex. 

All bridges considered so far had non-empty supports; however, this will not be the case in a 
typical graph. As an example, consider the edge of tenacity 17 in Figure 16. Since it does not 
assign minlevels to either of its endpoints, it is a bridge and it clearly has empty support. DDFS 
will discover this right away since the bud* of both of its endpoints is a. 

At the end of search level i = {t — l)/2, i.e., once MAX is done processing all bridges of tenacity 
t, compute the following set 

B = {bud*(t;) I tenacity (v) = t}. 

Then, with the help of the next lemma, all blossoms of tenacity t can be identified. The proof of 
this lemma is straightforward and is omitted. 

Lemma 29 The set of bases of all non-empty blossoms of tenacity t is precisely B. For each 
b & B, the set St defined in Definition 16, for blossom Bh^t is precisely {v \ tenacity(f) = 
t and bud*(w) = b}. 




Figure 17: Constructing a minimum length augmenting path between unmatched vertices /i 
and /2. 
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7.3 Finding the augmenting paths 



As stated at the beginning of this Section, during search level jm, where Im = ^jm + ^ is the length 
of a minimum length augmenting path in G, a maximal set of such paths is found; observe that 
Im is also the minimum oddlevel of an unmatched vertex. In contrast to previous search levels, 
in which all DDFSs will end with a bottleneck, in search level jm DDFS performed on certain 
bridges will end in two unmatched vertices. In this section, we give the procedure that is followed 
each time such a bridge is encountered; in effect, the operations of MAX have to be enhanced 
during search level jm- 

Besides the above-stated event, i.e., at a certain search level, DDFS ends in two unmatched 
vertices, another possibility is that G contains two or more unmatched vertices but it has no 
augmenting paths, i.e., the current matching is maximum (if so, no unmatched vertex will have 
a finite oddlevel). The algorithm will recognize this when it is done assigning minlevels and 
maxlevels to as many vertices as it could, i.e., it has explored the unmatched edges incident at all 
vertices having a finite evenlevel, the matched edge incident at all vertices having a finite oddlevel, 
and has performed DDFS on all bridges of finite tenacity that it has identified. 

7.3.1 Finding one augmenting path 

In Figure 15, when DDFS performed on bridge (u, v) of tenacity 11 at search level 5, it ends with 
the two unmatched vertices. The next task is to find an augmenting path of length 11 between 
the two unmatched vertices and containing edge {u,v). We will describe this via the graph in 
Figure 17. In this figure, minlevel(tt) > minlevel(t') and hence edge (u, v) is a bridge. DDFS on 
this bridge terminates with the two unmatched vertices, /i and f2- At this point, the stack of the 
DFS starting at u (v) contains all the vertices of tenacity 1^, that lie on the part of the augmenting 
path from u to fi (v to /2); observe that the bottom of the latter stack will contain b = bud(t;). 
All vertices of tenacity less than Im that constitute such an augmenting path are missing; they 
lie in petals whose bud*s, which are of tenacity Im, sit on the two stacks. The procedure given 
below will find one such complete augmenting path by recursively "opening" the nested petals. 

Let us show how to construct the path from v to f2- Since bud(i;) = b, we first need to find an 
evenlevel(6; v) path from v to b. The actions arc different depending on whether v is outer or 
inner; in this case v is inner. Therefore, evenlevel(6; w) = maxlevcl(6; t;), i.e., the path must use 
the bridge of this petal, which is (c, d). By jumping from v to its petal-node, the algorithm can 
get to the endpoints of this bridge. The "red" and "green" colors on the vertices of this petal, as 
assigned by DDFS (see Section 6), indicate that v was found via the DFS starting at c, say this is 
the red tree. The algorithm does a DFS on the red edges of this petal, starting from c, and finds 
a red path to v. Also, it does a DFS from d on the green edges to find a green path to b. 

By the rules set above, the DFS that started at d must have skipped to a = bud(tt;) while searching 
along edge {d,w). Therefore, the green tree yields the "path" d,a,b. At this point, we need to 
recursively "open" the petal whose bud is a and find an evenlevel(a; w) path in it. Again, we 
ask whether w is outer or inner. This time the answer is "outer" and so we simply keep picking 
predecessors of vertices until we get from w to a. This path is inserted in the right place in the 
"path" from d to b. The path from b to /2 is obtained via the same process: following down 
predecessors and recursively finding paths through any petals that are encountered on the way. 
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7.3.2 Finding a meLximal set of paths 

Section 7.3.1 gave the procedure for finding one minimum length augmenting path, say p. We 
now build on it to find a maximal set of disjoint minimum length augmenting paths; recall that 
this was the objective of a phase. We first show how to identify the set of vertices that cannot be 
part of an augmenting path that is disjoint from p. These will be removed from the graph, and 
MAX will proceed until it encounters another bridge which makes DDFS end with two unmatched 
vertices. The processes are repeated until all bridges of tenacity Im are processed. Since we do 
not remove any useful vertices, maximality is guaranteed. 

First, all vertices of p are removed; together with a vertex, all its incident edges are also removed. 
As a result, we may create vertices that have no more predecessors. Each such vertex v is also 
removed since the current graph does not contain a minlevel(i;) path anymore. This process is 
continued until every remaining vertex, other than unmatched ones, has a predecessor. 

Let us argue that the next time DDFS ends with two unmatched vertices, we are guaranteed 
to find an augmenting path between them. The main question is, "How are we sure that the 
procedure of Section 7.3.1 will be able to find appropriate paths through blossoms that have lost 
some of their vertices?" 

The structural properties established in Section 4 render the answer to this question surprisingly 
simple. Assume that vertex w G p is in blossom B and that vertex u is in the left-over part of this 
blossom. By Corollary 13, p must contain the iterated bases of v and the next minimum length 
augmenting path, using u, must also contain the iterated bases of u. However, the base, b, of 
blossom B is an iterated base of both the vertices! Therefore, b will be removed when p is removed 
and u cannot be on any minimum length augmenting path that is disjoint from p. Indeed, it is 
easy to see that the process of iteratively removing vertices that have no predecessors left will end 
up removing all of B. 

8 Proof of correctness and running time 

Lemma 30 Let {u,v) be a bridge such that tenacity (t;) = tenacity (u, f) < Im- Then the expres- 
sion for tenacity('u, v) uses minlevel(t;). 

Proof : If (u, v) is matched, u and v must both be inner vertices, since otherwise {u, v) will be 
a prop. Also, in this case, the expression for tenacity (tt, v) uses oddlevel(v), which is minlevel(v). 

Next assume that {u, v) is unmatched. Therefore the equality tenacity('u) = tenacity('U, v) implies 

oddlevel(v) + evenlevel(t') = evenlevel('u) + evenlevel('u) + 1 =^ oddlevel(n) = e{v) + 1. 

Now, if oddlevel(w) = minlevel(f ) then u must be a predecessor of v, contradicting the hypothesis 
that {u,v) is a bridge. Therefore, oddlevel('u) = maxlevel('y) . On the other hand, the expression 
for tenacity (-u, uses evenlevel('i;) and hence it uses minlevel(i;). □ 

Lemma 31 For each bridge (n, v) having tenacity 2i + 1 < Im, the algorithm will determine that 
{u, v) is a bridge and that its tenacity is 2i + 1 by the end of execution of procedure MIN at search 
level i. 
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Proof : First assume that (u, v) is matched. Since it is not a prop, u and v must both be inner 
vertices with oddlevel(?/) = oddlevel(u) = i. Therefore, during search level i, MIN will determine 
that {u, v) is a bridge and that its tenacity is 2i + l. 

Next assume that {u, v) is unmatched and that evenlevel(i() < evenlevel(f ) and that {u, v) is first 
scanned from u. Clearly evenlevel(M) < i. Since {u, v) is not a prop, v must already have a 
minlevel and minlevel(t;) < i. By Lemma 3, tenacity (t;) < tenacity(n, i;). Now if tenacity (t;) < 
tenacity(ti, v), both levels of v are already known and hence tenacity(u, f ) can be determined. If 
tenacity(t') = tenacity(ti, v) then by Lemma 30 the expression for tenacity(u, v) uses minlevel(i;) 
and hence tenacity('U, v) can be determined. Further notice that since tenacity (u, v) = 2i + 1, 
minlevel(i;) must be i. □ 

Theorem 32 For each vertex v such that tenacity(f ) < Im, Algorithm 1 assigns minlevel(i') and 
maxlevel(?;) correctly. 

Proof : We will show by strong induction on i that at search level f. Algorithm 1 assigns a 
minlevel of i + 1 to exactly the set of vertices having this minlevel and it assigns appropriate 
maxlevels to exactly the set of vertices having tenacity 2i + 1. The basis, i.e, i = 0, is obvious. 

To prove the induction step, assume that the hypothesis is true for all search levels less than i. Let 
minlevel(i;) = z + 1, let p be a minlevel(u) path, starting at unmatched vertex /, and let {u, v) be 
the last edge on p. It is easy to see that u must be BFS-honest w.r.t. p; if not, then v must occur 
on the shorter path to u, which contradicts the assumption that minlevel(u) = i + 1. Therefore, 
if \p[f to u]\ is odd (even), oddlevel(?x) = i (evenlevel(w) = i). By the induction hypothesis, u 
must be assigned this level, regardless of whether it is the minlevel or maxlevel of u. Therefore, 
on searching along an appropriate parity edge incident at u, MIN will find v. By the induction 
hypothesis, at the start of search level i, the minlevel of v is not set. Therefore, when edge (u, v) 
is examined, it is either still not set or it is set to i + 1 through some other edge scanned in the 
current search level (the latter case happens only if oddlevel(t') = i + Hence, while searching 
along edge (u, v), v will be assigned its correct minlevel, u will be declared a predecessor of v and 
(it, v) will be declared a prop. 

Next, assume t{v) = 2i+l and let p be a maxlevel(f ) path. By Theorem 25 there exists a unique 
bridge of tenacity t on p, say it is (a, h). By Lemma 31, by the end of MIN in search level i, (a, h) 
must be inserted in the list B{2i + 1). Therefore, MAX will call DDFS with this bridge to find 
its support. By the induction hypothesis, maxlevel(f) was not set in any of the previous search 
levels^. However, it may be set in the current search level; if so, v already belongs to a petal. If 
it is not set, DDFS will reach v and assign it its maxlevel. □ 

In Figure 18, the algorithm determines that {u, v) is a bridge of tenacity 15 at search level 6. 
However, according to Algorithm 1, DDFS has to be performed on {u,v) at search level 7. The 
question arises, "Why wait till search level 7; why not perform DDFS on (it, v) when procedure 
MAX is run at search level 6?" To clarify this, let us change the algorithm so it runs DDFS on 
an edge as soon as its tenacity and its status as a bridge have been determined. Assume further 
that among the various bridges ready for processing, ties are broken arbitrarily^. 

''The importance of this subtle point, which is related to the idea of "precise synchronization of events" is 
explained below with the help of Figures 18 and 19. 

*By making the example given in Figure 19 slightly bigger, one can easily ensure that there are no ties of this 
sort. 
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Figure 18: Can DBFS be performed on bridge Figure 19: If so, vertices a and b will get wrong 
{u, v) at searcii level 6? tenacities. 

Now consider the enhanced graph of Figure 19, in which vertices a and b are clearly in the support 
of the bridge of tenacity 13 and hence have tenacity 13. Assume that when MAX is run at search 
level 6 the bridge of tenacity 15 is processed first. Since the tenacities of vertices a and b are 
not set yet, DBFS will visit them and assign them a of tenacity 15, which would be incorrect. 
Observe that the correctness of MAX crucially depends on assigning tenacities to each edge of 
tenacity less than 2i + 1 before processing bridges of tenacity 2i + 1, i.e., the precise manner in 
which events are synchronized in Algorithm 1. 

Theorem 33 The MV algorithm finds a maximum matching in general graphs in time 0{m^/n) 
on the RAM model and 0{m^/n■ a{m, n)) on the pointer model, where a is the inverse Ackerman 
function. 

Proof : Through arguments made so far, it should be clear that each of the procedures of MIN, 
MAX, finding augmenting paths, and removing vertices after each augmentation will examine 
each edge a constant number of times in each phase. The only operation that remains is that of 
computing bud* during BBFS. This can either be implemented on the pointer model by using 
Tarjan's set union algorithm [Tar 75], which will take 0{m ■ a{m,n)) time per phase, or on the 
RAM model by using Gabow and Tarjan's linear time algorithm for a special case of set union 
[GT85], which will take 0{m) time per phase. Since 0{^/n) phases suffice for finding a maximum 
matching [HK73, Kar73], the theorem follows. □ 

Theorem 33 raises the question of finding a linear time implementation of bud* on the pointer 
model; this was left as an open problem in [Vaz94]. One way of accomplishing this is to prove the 
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claim made in [MV80] that because of the special structure of blossoms, path compression itself 
suffices, together with a charging argument that assigns a constant cost to each edge. 

9 Equivalence of definitions 

Below we establish equivalence of the two definitions of blossoms. Let us start by providing the 
definition of blossom as given in [Vaz94] ; we will denote such a blossom of tenacity t < Im and 

base b by Bf,^t"- Let f be a vertex with tenacity(t;) < t. We will say that an outer vertex b is 
base>t(u) if for some positive k, base'^(f) = b, tenacity(6) > t, and tenacity(base'^~^(?;)) < t. 
Then 

Sb,t° = {v I tenacity (-y) < t and base>t('y) = b}. 
Proposition 34 The two definitions of blossom are equivalent, i.e., Bb,t = Bb,t"- 

Proof : Let v € Bh^t- We will show that v G Bb,t° by considering the following three cases. The 
set St is defined in Definition 16. 

1. V e St. In this case, tenacity(t;) = t and base(t;) = b, and therefore base>t(i;)) = b. Hence 

V e Bb,t°. 

2. t> G Bb^t-2- In this case, tenacity (w) < t and for some k > \, base*^(i;) = b. Clearly, 
base^~^(i') G Bb,t-2 and therefore base>t(i;) = b. Hence v G Bb,t° ■ 

3. G Bu,t-2 and u G St. In this case, tenacity(w) < t, tenacity(u) = t, base('u) = b, and for 
some k >1, base^(i;) = u. Therefore, hase'^~^^{v) = b and base>t(u)) = b. Hence v G Bb,t°- 

Next, let V G Bb^t°- Once again we will consider three cases to show that v G Bb^f 

1. tenacity(u) = t. In this case, base(f ) = b and therefore v G Sf. Hence v G Bb,t- 

2. tenacity('u) < t, for some k > 1, base'^(t') = b and tenacity(base'^~^(v)) < t. In this case, 

V e ^6,4-2. Hence v G Bb,t- 

3. tenacity(u) < t, for some k > 1, base'^(f ) = b and tenacity(base'^~^(u)) = t. Let base^~^(f ) = 
u. Then, tenacity('u) = t and base('u) = b. Therefore, u e St and v G Bu,t-2- Hence v G Bb,t- 

□ 
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